First excitations of the spin 1/2 Heisenberg antiferromagnet on the kagome lattice 
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We study the exact low energy spectra of the spin 1/2 
Heisenberg antiferromagnet on small samples of the kagome 
lattice of up to TV = 36 sites. In agreement with the conclu- 
sions of previous authors, we find that these low energy spec- 
tra contradict the hypothesis of Neel type long range order. 
Certainly, the ground state of this system is a spin liquid, but 
its properties are rather unusual. The magnetic (AS = 1) 
excitations are separated from the ground state by a gap. 
However, this gap is filled with nonmagnetic (AS = 0) exci- 
tations. In the thermodynamic limit the spectrum of these 
nonmagnetic excitations will presumably develop into a gap- 
less continuum adjacent to the ground state. Surprisingly, 
the eigenstates of samples with an odd number of sites, i.e. 
samples with an unsaturated spin, exhibit symmetries which 
could support long range chiral order. We do not know if these 
states will be true thermodynamic states or only metastable 
ones. In any case, the low energy properties of the spin 1/2 
Heisenberg antiferromagnet on the kagome lattice clearly dis- 
tinguish this system from either a short range RVB spin liquid 
or a standard chiral spin liquid. Presumably they are facets 
of a generically new state of frustrated two-dimensional quan- 
tum antiferromagnets. 



I. INTRODUCTION 

Zeng and Elser jij] were the first to point out that the spin 
1/2 Heisenberg antiferromagnet on a kagome lattice (hereafter 
called KHA) could be disordered. Since then many studies 
have been devoted to this system. All quantum approaches 
based on either exact diagonalizations j^J^| , on perturbational 
series expansions [^] or on high temperature expansions 
point to a disordered ground state. Previous exact diago- 
nalization studies have shown that spin-spin, spin ne- 
matic, spin Peierls and chiral-chiral correlations decrease very 
rapidly with distance. Therefore the idea that the KHA is a 
spin liquid is by now widely accepted. The precise nature 
of the spin liquid state of the KHA and of its low lying ex- 
citations is, however far from clear. Some approaches point 
to a liquid of short range dimers |]-^|, others suggest that 
the KHA might be a candidate for a chiral spin liquid state 
of the kind first proposed for the triangular antiferromagnet 
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by Kalmeyer and Laughlin P |ll"l. It seems to be the com- 
mon view that in both pictures any excitations, magnetic or 
nonmagnetic ones, are separated from the ground state by a 
finite gap. While for the magnetic excitations this view is sup- 
ported by all previous numerical work on the KHA Jl|Q, the 
situation is less clear for the nonmagnetic excitations. The 
only piece of evidence for a finite albeit small gap separat- 
ing the singlet excitations of the KHA from its ground state 
comes from the work of Zeng and Elser ||. Elaborating on 
the idea of a short-range dimer liquid these authors construct 
an effective low energy hamiltonian for the singlet subspace 
of the KHA, which as expected, exhibits a small gap above 
its ground state. The question of whether one of these pic- 
tures, that of a short range RVB spin liquid or that of a chiral 
spin liquid, applies to the KHA has led us to study not only 
the energies but also the symmetry properties of a very large 
number of low lying levels of the exact spectra of the KHA. 



II. NUMERICAL APPROACH AND SOME 
GENERAL RESULTS 

Using the Lanczos technique and a complete group theoret- 
ical analysis (as in we computed the low energy spectra 
of the Heisenberg hamiltonian 



T~L = J ^ 2Si.Sj 



(1) 



<ij> 



for small samples of the kagome lattice on a 2-dimensional 
torus with periodic and twisted boundary conditions (in 
the following J = 1). We have obtained a very large 
number of low lying levels in each irreducible representa- 
tion (IR) of SU{2) and of the space groups of the N = 
9, 12, 15, 18, 21, 24, 27, 36 samples. 

As in all spectra of Heisenberg antiferromagnets that we 
have examined, the ground state energy of the S-subspaces 
increases with the total spin S. However, the spectra of the 
finite samples of the KAH do not show the pattern that is 
characteristic for systems which exhibit Neel type long range 
order in the thermodynamic limit as, e.g. the triangular anti- 
ferromagnet (for a detailed discussion of the signature of long 
range order in the spectra of finite samples see ref. |l^l). In- 
stead, all the signatures of a "liquid" are present in the low en- 
ergy spectra of the KAH: i) The lowest levels associated with 
the different momenta k in the Brillouin zones of the samples 
are almost independent of k. For instance, for the largest 
sample (TV = 36), the energies of these lowest levels vary by 
less than 0.2% when k is varied through the four inequivalent 
k-points of the Brillouin zone (see Table 1). This absence of 
dispersion excludes the possibility of a broken translational 
symmetry, ii) The system is extremely soft against any twist 
of the boundary conditions . Hi) The size dependence of 
the ground state energy is more than one order of magnitude 
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smaller than for the triangular antiferromagnet in the same 
range 9 < N < 36. 



III. MAGNETIC GAP AND NONMAGNETIC 
EXCITATIONS 

Our results, (Fig. 1), confirm the conclusion of Elser and 
coworkers [jj],^] about the existence of a gap for the magnetic 
AS — 1 excitations. The thermodynamic limit of this gap 
cannot be extracted from the present data with high accu- 
racy. The extrapolations shown in Fig. 1 point to a lower 
bound of about J/20 (see also the discussion in ref. In 
any case, the value of the magnetic gap appears to be at least 
one order of magnitude smaller than the exchange energy J 
which is needed to break an isolated singlet pair. This is a 
first indication that the picture of the kagome antiferromagnet 
as a spin liquid consisting of short-range singlet dimer pairs 
may be inappropriate. Nevertheless, it comes as a surprise 
that quite contrary to the standard picture of a spin liquid 
the magnetic gap of the KHA is filled with nonmagnetic ex- 
citations. This is illustrated in Fig. 2, where we display the 
integrated density of states of the N — 36 sample. In this 
sample there are 183 singlet levels below the lowest triplet. 
Similarly, in the N = 27 the AS = 1 gap is filled with 153 
nonmagnetic excitations. For the samples we have examined, 
the number of nonmagnetic (AS = 0) excitations within the 
magnetic gap grows roughly as a N with the system size, where 
a ~ 1.15 and a ~ 1.18 for the even and odd samples, respec- 
tively. These results call for a few comments: 

i) The density of nonmagnetic levels above the ground 
states of the two largest samples (N = 27, 36) strongly sug- 
gests that in the thermodynamic limit the nonmagnetic exci- 
tation spectrum of the KAH is a gapless continuum adjacent 
to the ground state. Interestingly, an algebraic decay of cer- 
tain correlation functions would be compatible with such a 
gapless nonmagnetic continuum. A candidate is the dimer- 
dimer correlation function, 

C<!,]Xt,i> =< (Si • Sj)(Sfe • Si) > — < Si • Sj > 2 , (2) 

which might decrease algebraically with the distance between 
the nearest neighbor bonds < i,j > and < k, I >. Leung and 
Elser ^ have calculated this function for the ground state of 
the N = 36 sample of the KHA. Their result, a damped os- 
cillatory decrease of C < ij >< k,i>, is not inconsistent with an 
algebraic decay. An algebraic decay of the dimer-dimer corre- 
lations would imply that the system is critical with respect to 
nonmagnetic quantum fluctuations. Whether this is the case 
for the KHA or not will be impossible to decide from finite 
size studies of the correlation function alone. 

ii) The exponentially large number of low lying singlets 
of the spin 1/2 KHA is reminiscent of the ground state de- 
generacy of the corresponding classical model. Stability con- 
siderations show that only planar spin configurations qualify 
as true classical ground states ]7|Jl^-p^]. Their number and 
hence the ground state degeneracy of the classical KHA grows 
as 1.134^ with the system size N & The different planar 
configurations are connected with each other by local rota- 
tions in spin space Ufx- In a semiclassical picture these local 



rotations provide tunneling paths between the different pla- 
nar configurations. It is therefore tempting to think of the low 
lying singlets of the quantum KHA as of tunnel-split classical 
ground states. This semiclassical picture has been pursued 
by von Delft and Henley |1S[ ]. As one of their main results 
these authors find that for spin 1/2 all the tunneling events 
involving the coherent rotation of the spins on simple loops 
on the kagome lattice yield zero tunneling amplitudes due to 
destructive interferences between different tunneling paths. 
The smallest effective tunneling event requires the coherent 
rotation of the spins on two nested hexagons on the kagome 
lattice involving 24 spins in total. Certainly, this scenario 
cannot explain the abundance of low lying singlets in a sam- 
ple consisting of only 21, 27, or 36 spins. In our view, this 
high density of low lying singlets in small samples of the KHA 
is a strong argument against the validity of the semiclassical 
approach for the spin 1/2 KAH. 

Hi) Another comparison that suggests itself is between the 
density of low lying singlets and the dimension of the near- 
est neighbor valence-bond basis on the kagome lattice which 
grows as 2 JV/3 = 1.26^ with the system size §. From this 
point of view the variational approach of Zeng and Elser |]] 
which builds on the valence bond basis seems fully justified. 
The main difference between the approximate singlet spec- 
trum of Zeng and Elser || and our numerically exact result 
lies in the density of levels at the bottom of this spectrum 
(see Fig. 2). The absence of a gap above the ground state 
and the high density of very low lying levels in the exact 
spectrum suggests that a dimer product representation of the 
corresponding eigenstates will necessarily contain long range 
singlet pairs. The importance of such longer range singlets 
in the spin liquid picture of the KHA might have been an- 
ticipated from the work of Zeng and Elser who observed that 
the inclusion of first and second-neighbor singlet pairs in their 
variational Hilbert space led to a considerable improvement 
over the results obtained in a pure first-neighbor dimer basis. 

iv) The exponential number of non magnetic states in the 
gap should be visible in various experimental situations: It 
may explain the vanishing of the neutron elastic forward scat- 
tering cross section |2(J an d the very weak field dependence 
of the low temperature heat capacity of SrCrGaO observed 
by Ramirez and coll. ^l| . Specific heat data and neutron 
scattering experiments both point to a density of states that 
increases linearly with the energy. In the very low energy re- 
gion (E/N = —0.4384... — 0.4364) our numerical data are 
consistent with such a linear energy dependence of the den- 
sity of states (see the inset in Fig. 2, where the integrated 
density of states is plotted); they are definitely not consistent 
with an energy independent density of states as it would be 
obtained for an ordinary spin glass p2| . 

IV. SPIN 1/2 EXCITATIONS 

A. Symmetries of the low lying S=l/2 eigenstates. 

The second surprise lies in the degeneracies of the low ly- 
ing levels of the spectra of the samples with an odd number 
of spins (N — 9, 21, 27) and in the symmetries of the eigen- 
states associated with them. We are focusing here on the ex- 
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ponentially large number of spin doublets which fill the gap 
between the ground state and the lowest 5 — 3/2 eigenstate. 
On account of the finite size scaling of the available spectra, 
we can ascertain that in the thermodynamic limit all these 
eigenstates are complex and exhibit an at least twofold de- 
generacy in addition to their magnetic degeneracy. (In the 
following we shall consistently ignore the trivial spin degener- 
acy.) Hence all these states may participate in a spontaneous 
mechanism breaking the time reversal (and parity) invariance 
of the Heisenberg hamiltonian. In the spectra of the smallest 
samples, there are still some low lying levels which belong to 
ID IRs of the space group and are hence nondegenerate, but 
when the size of the sample increases these levels are pushed 
towards higher energies: in the N — 27 sample there are 153 
5=1/2 levels below the lowest 5 = 3/2 level, only 3 of them 
belong to ID IRs of the space group. The finite size extrapola- 
tion predicts that for odd sizes N > 50 no states from the ID 
IRs of the space group will be left in the magnetic gap. This 
is indeed a very specific property of the KAH which we have 
never observed in the spectrum of any spin hamiltonian on 
the triangular lattice : in the usual triangular Neel antiferro- 
magnet the homogeneous low lying states belong to the trivial 
representation of the Cz group, whereas in the present case 
the 5=1/2 levels of the trivial representations are pushed 
high up in the spectra.^] This observation leads us to con- 
centrate on possible new properties of this low lying 5=1/2 
continuum. 



B. Chiral observable and correlation function 

All these low lying levels (except the 3 belonging to ID 
IRs) can sustain non-zero expectation values of the chirality 
operator 

E 123 = 2 Si ■ (Sa x S 3 ) = Im(P 123 ) (3) 

where P123 is the cyclic permutation of three spins around a 
triangular plaquette. 

We have measured this quantity in the lowest lying states 
of the N = 9, 21 and 27 samples for triangles with side lengths 
1, \/3, 2, y/7, 3 (see Table 2). Compared with the eigenvalues 
±V3/2 of the operator E123 in the 5=1/2 bound states 
of three spins, the measured expectation values are indeed 
quite small, and there is a sharp decrease from the N = 9 
to the N = 21 sample. But in view of the data for N = 21 
and -/V = 27 it is unclear how this observable will behave for 
larger sizes. The chiral-chiral correlation function ( last lines 
of Table 2) is indeed short ranged, which may explain the 
quasi-absence of size effects between N = 21 and iV = 27, 
and is indeed characteristic of a liquid. 



1 In the case of the q = Neel order, enforced on the kagome 
lattice by a large enough second neighbor antiferromagnetic 
coupling |^3| , the ground state of the odd samples does in fact 
belong to the two complex conjugate IRs of C3 [^3). However 
in this case it is the only low lying state with this symmetry, 
all the others have much higher energies. 



In order to have an indication of the possible order of mag- 
nitude of the chirality in a pure chiral spin liquid, we have 
computed the expectation value of E123 for the spin 1/2 exci- 
tations of the Laughlin wave function jll| on the same small 
lattices; these expectation values are denoted (L.w-f) in Table 
2. While for triangles of size 1 the L.w-f values are an order 
of magnitude larger than the exact results, they are of the 
same order of magnitude for larger sizes. By this comparison 
we do not want to suggest that the physics of the KAH is the 
same as that of the Laughlin picture of a chiral spin liquid ( 
as has been found in section 4. above, the first excitations are 
different). We only want to point out that the exact results 
may be significant in spite of their numerically small values. 

On the basis of results obtained by applying the 517(2) 
Schwinger- boson approximation to samples of the KAH with 
sizes of up to N = 72 sites Sachdev argues that the chirality 
of the ground state goes to zero with increasing system size 
It is indeed true that the exact ground states of the even 
N samples have zero chirality for N = 12, 18, 24, 36. The 
present work indicates that this may not be true for the odd 
N samples. 



C. Chern index of the low lying S=l/2 levels 

Searching for a quantity that might provide further ev- 
idence for or against collective chiral behavior we followed 
the proposal of Haldane and Arovas |m| and computed the 
Chern number of the homogeneous (k = 0) ground state 
of the N = 9 sample. It is a state which transforms as 
Rz-k/s^ — e %2lT ^ J il) under a spatial 27r/3 rotation. In this en- 
tire paragraph we shall concentrate on these specific states, 
which are the low lying homogeneous states of the 5=1/2 
continuum that we have been discussing above. Anticipating 
further results we shall refer to these eigenstates as "chiral" 
states henceforth. In the N = 27 sample, 36 of the 39 homo- 
geneous states of this low lying continuum are chiral states, 
and we expect that all states of this continuum will be chiral 
in the thermodynamic limit. The Chern number is a topolog- 
ical property of the eigenfunctions of a system constrained by 
twisted boundary conditions J24|. If Ti and T2 denote the 
two vectors defining the sample cell, the twisted boundary 
conditions are defined by 

5 ± (R, + T J ) = e ±i *'S ± (R 4 ), (4) 

where $1 and $2 are the two angles defining the rotations 
of the spins around the z direction. Let us denote by 
|n($i, $2) > the eigenstates obtained by adiabatic evolution 
of the state \n > under the twists. The Chern index of the 
ket \n > is defined by 

It is a topological constant that can take only integer values. 
It can only be non-zero for complex, degenerate eigenstates, 
and, as will be discussed below, a non-zero value of the Chern 
number of an eigenstate \n > implies peculiar physical prop- 
erties of the system in this state. A formula equivalent to (4) 
is obtained by the use of Stokes' theorem: 



3 



C(n) = 2ni <f> < n\dn > 



(6) 



It shows that the Chern number counts the number and the 
nature of the singularities of the phase of |n($i, $2) > in the 
($1, $2) Brillouin zone. It is a measure of the vorticity of 
this phase. For computational purposes the expression (4) is 
quite cumbersome. A more convenient expression is obtained 
by changing to a spatially varying reference frame for the spin 
variables such that the twisted boundary conditions (3) are 
replaced by periodic ones. The main steps in this procedure 
are as follows: Let So be the reference frame for the spin at 
the origin Ro of the sample. Then the reference frame at the 
site Ri is obtained from So by a rotation through the angle 



B(Ri) = (R - R )(ei6>i + e 2 6 2 ), 



(7) 



around the z-direction, where ei, e 2 span the Bravais lattice 
and 61, 62 are the increments of the rotation angle along ei 
and e2. Oi, 62 are chosen such that they add up to the twist 
angles, i.e.: 

$ a = 6(R + T a ), a = 1,2. (8) 
In this spatially varying reference frame the hamiltonian reads 

H{e u 0a) = 2 J2 i 3 * 1 ^ + V 2 (e x -iS+S~ + h.c.)}, (9) 



where 



<i,3> 



0(Ri) - 0(Ri 



(10) 



The Chern number in the eigenlevel \n > associated with the 
eigenvalue E n of TL is then readily obtained as the average 
over the ($1, $2) Brillouin zone of the function : 



K a t, = 2ne ab 



(E n — E p ) 



where 



h.c. 



(11) 



(12) 



The phase angle Xij depends linearly on 3> a through the lo- 
cal twists (7), and is a linear combination of the space 
components of (R, — Rj). 

The computation of the Chern number is a heavy task. 
After a thorough study of the chiral ground state of the N = 9 
sample we can ascertain that for this state the Chern number 
is +1 (respectively —1 for the complex conjugate eigenstate). 
In view of our study of the singular points of E„(&i, $2) in 
the chiral states of the N = 21 and 27 samples, we expect it to 
be odd for any of the low lying doublets and most probably 
equal to ±1 for the ground state (in this last case, we find 
only one conical point at ($i,$2) = (0,0)). 

The physical significance of the Chern index may be in- 
ferred from the expression (11) which is a zero frequency Kubo 
response function. The general formula for the response of an 
observable B to an excitation of the observable A reads at 
T = 0: 



Xba(<*> = 0) = 2h ImJ2 



< n\B\p >< p\A\n > 



pi^n 



{En 



Epf 



(13) 



where A is the time derivative of the observable coupled to the 
external field. Comparing this general form of the response 
function with the special form (11), one sees that in (11) 
both operators A and B are linearly connected with the total 
spin currents J a = *X/<ij> { eXij ^t^j~ — h.c.), where the 
^2<ij> * s ^° ^ e restricted for a = 1,2, to the bonds along 
the directions ei and &2- The external field that drives a spin 
current along, e.g. the ei direction, is a magnetic field with a 
constant gradient in this same direction. This is most easily 
seen on the example of the square lattice, where the relation 
between (61,62) and ($1,^2) is diagonal. The case of the 
kagome lattice is technically slightly more involved, but the 
physics is the same. The perturbation of the hamiltonian H 
induced by a magnetic field B z with a constant gradient in 
the ei direction is given by 



V = 



ctiXuSi , 



(14) 



where Xn = (Rj • ei). Here, ct\ is the linear increment of 
the magnetic field over one lattice constant. (Here and in the 
sequel we set ft, the Bohr magneton and the gyromagnetic 
factor equal to unity). Identifying the observable A as A = 
J^(Ri • ei)^ 2 one sees immediately that its time derivative is 
the total spin current in the ei direction: 

A = -i [A, H] = Im{e zx ^SfSj}(X rl - X 3l ) = Ji 



<i,3> 



(15) 



In the orthogonal geometry one has the simple relation J a = 
■§§-, a = 1,2, and the expression (11) takes the form of a 
transverse current- current correlation function: 



2 Im 



\ - < raj J a \p >< p\J b \n > 



p^n 



{E n 



E P ) 2 



(16) 



This correlation function measures the transverse current J2 
generated in the e2 direction by a gradient of the magnetic 
field B z along the ei direction. Hence we arrive at the linear 
relation 



< J 2 >= 27rK a , 



dB z 

dxl 



(17) 



At this point it is important to note that the above derivation 
of (17) is purely formal. It is a linear expansion with respect 
to the perturbation V, (14), which, as it is obvious from (14), 
becomes arbitrarily large in the thermodynamic limit. Thus, 
(17) cannot be considered a physically meaningful relation be- 
tween the spin current and the gradient of the inhomogeneous 
magnetic field. However, by employing a gauge argument Hal- 
dane and Arovas |^j show that a physically sound relation is 
obtained from (17), if the response function K is replaced by 
its average over the torus < $1 , $2 < 27r of the twist angles, 
i.e. by the Chern number C(0). 

From the mathematical point of view the Chern number 
of this spin system is equivalent to the TKNN index of the 
Quantum Hall Effect (QHE) M, and as in the QHE the spin 
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1/2 excitations are separated from the S = ground state by 
a finite gap (see Ag =1 / 2 in Fig. 3). But there are very definite 
and crucial differences: i) In our spin system, contrary to the 
QHE, the parity and time reversal invariance are not exter- 
nally broken. Eigenstates with positive and negative Chern 
number are degenerate, ii) In each spin sector there is a con- 
tinuum of excitations adjacent to the ground state and there 
may be couplings between the two sectors under the effect of 
an external magnetic field. So it is difficult to imagine how the 
microscopic rigidity associated with the Chern number could 
become manifest on a macroscopic level, Hi) Finally our nu- 
merical results seem to indicate that the creation of two spin 
1/2 excitations is less favorable than the creation of a spin 1 
excitation : 2 •Eas=i/2 > Eas=i (see Fig. 3), but because 
of the uncertainties in the extrapolation procedure, this last 
result is not entirely reliable. The question of whether the 
excitation of two spin 1/2 entities is energetically more or less 
favorable than the creation of one S = 1 excitation is indeed 
an major open point. In the first case the spin 1/2 excita- 
tions would be true thermodynamic excitations, in the oppo- 
site case they could only appear as metastable excitations in 
sophisticated dynamic experiments. Whatever the answer to 
this last question may be, the KAH is certainly not a system 
which fits easily into the frame of standard continuum chiral 
theories. In any case the symmetries of these spin 1/2 exci- 
tations are features that will survive in the thermodynamic 
limit although the states may be metastable in this limit. We 
may hypothesize that their non-zero Chern number, since it 
is a quantum number, will also survive in the thermodynamic 
limit. In this picture the chiral S = 1/2 states of the KAH 
can certainly not be viewed as simple bound states of three 
spins 1/2 in a sea of singlets: a cluster of three spins in its 
S — 1/2 ground state has indeed a nonzero expectation value 
of the chirality H123 = ±\/3/2, but its Chern index is zero. 
The Chern number is a measure of a topological rigidity of the 
N-particle states, which is by itself unique. To our knowledge 
the KAH is the first system with a hamiltonian that does not 
break parity and time reversal invariance, where a non-zero 
Chern number has been observed. 



V. CONCLUSION 

In conclusion, our numerical study of the low lying spec- 
tra of the spin 1/2 KAH leads us to assert that this system 
is a strange spin liquid: Its lowest excitations are soft non- 
magnetic excitations whose existence in a spin system has not 
been anticipated in previous theoretical investigations of anti- 
ferromagnets. The absence of a gap between these excitations 
and the ground state may signal that the KAH is, at T = 0, in 
a critical state with arbitrarily long-ranged singlet- singlet cor- 
relations. In any case, the existence of this continuum of low 
energy nonmagnetic modes in the KAH already contradicts 
the conception of the ground state of this system as a short- 
range RVB state. It also contradicts the classical and semi- 
classical pictures of the KAH which are all more or less based 
on the idea of preferred planar spin arrangements, which im- 
ply a breaking of the SU (2) symmetry. As has been shown in 
previous work |l2[|l3] ] , the spectra of finite samples of a spin 
system have to exhibit a certain signature, if a breaking of 



the SU (2) symmetry is to occur in the thermodynamic limit. 
The spectra of the finite samples of the KAH do not show this 
signature [juj, and thus a planar arrangement of the spins of 
the KAH can be ruled out. 

In the sea of long-range correlated singlets an unpaired 
spin 1/2 is surrounded by transverse spin currents, which are 
another signature of long range correlations. The non zero 
Chern number of the 5=1/2 states of the KAH is a proof of 
some kind of microscopic rigidity (sensitivity to the bound- 
ary conditions), which is usually found in chiral theories in 
the literature. However, in standard chiral theories these spin 
1/2 excitations are indeed the lowest excitations of the sys- 
tem, which is definitely not the case for the KAH. We do 
not know if these spin 1/2 excitations will be true thermo- 
dynamic excitations or only metastable ones, but they might 
show up in dynamical experiments (transverse spin diffusion, 
spin echoes). 

To conclude we want to emphasize that there is some evi- 
dence that this kind of ground state and low lying excitations 
( critical singlets and gapped magnetic excitations) are robust 
properties of the KAH: finite perturbations (second neighbor 
interactions) are needed to drive the KAH into Neel ordered 
states Jl^,^], and the XY model on the kagome lattice has 
a spectrum that is similar to the spectrum of the Heisenberg 
model. Thus, the above properties of the KAH with nearest 
neighbor interactions may be facets of a new state of frus- 
trated quantum antiferromagnets which is generically differ- 
ent from the Neel ordered state on the one hand and from the 
short-range RVB state on the other |pr^-pl|. The low energy 
physics of the KAH may be embodied in real compounds: In 
spite of the fact that SrCrGaO is a spin 3/2 compound with 
many complexities (dilute triangular planes of spins interca- 
lated between the kagome planes, role of the defects), some 
of its unusual properties (low temperature specific heat, dy- 
namical spin correlations) may be explicable on the basis of 
the present results. 
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FIG. 1. a) Energy per bond j) in the ground states 

of S m i n and S m i n + 1. The energy per bond in the thermo- 
dynamic limit is approximately 2 < S^.Sj >= —0.43. b) Spin 
gap versus 1/N for the even TV and the odd N samples. 
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FIG. 3. Sas=i/2 : energy for creating one S =1/2 exci- 
tation. This energy is the difference between the interpolated 
ground state energies for S — 1/2 and S — (see Fig.l). 
Eas=i '■ energy for creating one S = 1 excitation in the even 
samples (magnetic gap). The graphs shown in this figure 
are based on polynomial interpolations between the numer- 
ical data. Continuations of these graphs to larger values of 
N (N > 27 for E AS=1/2 and N > 36 for E AS =i) are unwar- 
ranted. 
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FIG. 2. Integrated density of states of the low lying levels 
of the N = 36 sample plotted versus the energy per bond 
(heavy line). The light and the broken line display the same 
quantity for the S — and S — 1 subspaces separately. The 
inset shows a quadratic fit to the lowest 70 states of the spec- 
trum. 



TABLE I. Lowest energies per bond associated with the 
four different k-points in the Brillouin zone of the JV = 36 
sample. 
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(0,0) 
(f ,0) 

(f ,0) 


-0.438377 
-0.437851 
-0.437585 
-0.438096 



G 



TABLE II. Expectation values of the chiral operator for 
triangular loops of various sizes: d is the side-length of the 
triangle in units of the nearest neighbor distance. These ex- 
pectation values have been computed in the chiral ground 
states of the odd samples ( they are homogeneous k = (0, 0) 
states and belong to the two degenerate complex conjugate 
IRs of C3). The lines with the entry (L.w-f) in the left col- 
umn contain the expectation values of the same operators in 
the variational Laughlin wave- function of ref. O]. The last 
six lines contain the correlations between the chiralities on 
elementary triangular plaquettes for various distances. 
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